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ABSTRACT 

We present the detection of sodium absorption in the atmosphere of the extrasolar 
planet WASP- 17b, an inflated 'hot- Jupiter' in a tight orbit around an F6 dwarf. In- 
transit observations of WASP- 17 made with the MIKE spectrograph on the 6.5-m 
Magellan Telescope were analysed for excess planetary atmospheric absorption in the 
sodium I 'D' doublet spectral region. Using the interstellar sodium absorption lines as 
reference, we detect an excess 0.58 ±0.13 per cent transit signal, with 4.5cr confidence, 
at 1.5 A bandwidth around the stellar sodium absorption feature. This result is con- 
sistent with the previous VLT detection of sodium in WASP- 17b, confirming that the 
planet has a highly inflated atmosphere. 
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1 INTRODUCTION 

Numerous techniques have now been developed to study 
the atmospheric properties of transiting planets. For ex- 
ample, emission from the planet inferred by observations 
of the secondary eclipse can reveal t heir day-side temper - 
atures (e.g. Charbonneau et al. 20^; IPeming et al.lliooih . 
Additional infrared photometry of the phase variation can 
map the temperature di stribution of 'hot- Jupiter' plan- 
ets ([Knutson et al.l l2007l ) . As light from the host star 
passes through the upper atmosphere of a planet dur- 
ing the primary transit, it is possible to investigate the 
planetary atmospheric compo s ition via transmission spec- 
troscopy. ICharbonneau et al.l (|2002l ) employed this tech- 
nique to detect the atmosphere of an extrasolar planet for 
the first time. Their observations of HD 209458b with the 
Hubble Space Telescope (HST) revealed a 0.02 per cent tran- 
sit signal in the Nal resonance lines, indicative of Nal ab- 
sorption in the upper atmosphere of the planet. 

The amount of absorption measured by transmission 
spectroscopy probes the wavelength dependence of the 
height of the planet's atmosphere. The optical depth is de- 
pendent on the temperature, surface gravity, and compo- 
sition of the atmosphere. For 'hot-Jupiters', the opacity is 
predicted to be strongest at the Nal and KI resonance lines 
in the optical, whilst molecular species such as water and 
methane dominate the near- infrared part of the spectrum 
(|Seager Sasselovll2000l : lBrownll200ll : iHubbard et al ] |200lh . 



This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile, 
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In recent years, the use of transmission spectroscopy have 
lead to a series of discoveries in the study of 'hot- Jupiter' 
atmospheres, including t he escaping atomic hydro gen enve- 
lope around HD 209458 (IVidal-Madiar et al." 2003), and the 
presence of haze in HD 189733b ^Ppnt et al. 2008). 

Despite the recent successes, the expected planetary 
absorption signal is very small, and the technique of 
transmission spectroscopy remains difficult. So far, Nal 
detection via ground based high-resolution spectroscopy 
has only been achieved on t hree transiting planets, HD 
209458b (ISnellen et al.l l2008l : iLandand-Shula et al.l l2009l : 
Jensen et al. 2 0JT|), HD 189733 b(| Redfield et al.ll2008h . and 
WASP-17b (W ood et al.ll201lh. alone: with numerous un- 
successful attempts fe.g. lNarita et al.ll2005l : iBozorgnia et al.l 
I2OO6I ). Recent detections around the KI resonance line by 
nar rowband photom etry has also been achieved with XO - 
2b ("Sing et al.l lioT? ) and HD 80606b (jColon et al.ll2012l ). 
High-resolution spectroscopic detections often require 8- 
10 m class telescopes, with data taken in optimal condi- 
tions, yielding signal-to-noise ratios S/N ^ 100 pixel" ^. 
Even so, the transit signal is heavily contaminated by much 
larger systematic tr ends, such as echelle blaze funct ion vari- 
ation s (jWinn et al.lt2004 ). CCD non- li nearity ( Snellen et al.l 
2008h, and telluric contamination (e.g. Albrecht et al.l 20091 : 



Langland-Shula et al.l 1200*91 : I Jensen et al.l 2011). Here, we 



present a novel approach of using the nearby interstellar ab- 
sorption lines to correct for the majority of these systematic 
effects. 

WASP-17b (jAnderson et al.l l2010h is a 0.49 Mjup, 
1.99 -Rjup planet in a retrograde, 3.7 day orbit ('Bavliss e t al I 
I2OIOI : iTriaud et al.l 1,2010). The equilibrium temperature of 
the planet, measured by Spitzer photometry of the sec- 
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ondary eclipse, is 1771=b35K, assuming zero albedo and full 
re distribution (lAnderson et al. 2011). Recent observations 
bv lWood et alJ (|201lh with GIRAFFE on the VLT detected 
a 1.46 per cent transit in the Nal 'D' lines at 0.75 A band- 
width, a level ten times larger than detections made around 
any other planet. 

.Bavliss et al] (|2010l ) obtained time-series spectra of 
WASP- 17 during the primary transit event to measure the 
spin-orbit alignment of the system. In this study, we present 
an analysis of this data to monitor for Nal absorption in the 
planetary atmosphere. 



2 OBSERVATIONS AND DATA ANALYSIS 

Observations of WASP- 17 were obtained with the Magellan 
Inamori Kyocera Echelle (MIKE) spectrograph on the 6.5- 
m Magellan II (Clay) telescope, located at Las Campanas 
Observatory, on the night of 2010 May 11. The spectra were 
obtained at a resolution of A/ A A ~ 48000 in the wavelength 
range 5000-9500 A, over echelle orders 37-69. Each exposure 
was 600 s in length, giving an average S/N of 82 pixel" ^ over 
the Nal 'D' doublet region, in order 58. Th-Ar arc lamp 
exposures were taken bracketing each exposure. Of the 37 
object exposures, 21 were taken out-of-transit, 11 in-transit, 
2 during ingress, and 3 during egress. Two observations were 
rejected for having S/N < 60, one was taken during egress, 
the other whilst out-of-transit near the end of the night, 
at high airmass. Observations were paused for 30 minutes, 
whilst the planet was in-transit, due to the zenith limit of the 
telescope. T he details of the obse rvations and data reduction 
are given m iBavhss et aP (|201Q| ). 



2.1 Spectral analysis 

2.1.1 Correction of systematic trends 

A wavelength shift of ^ lkms~^ was registered over 
the night, due to a combination of instrument drift, air- 
mass change, system orb ital dynamics , and the Rossiter - 
McLaughlin (RM) effect ( Ros siterlll924l : "McLaughlin ' 1924') . 
The average shift per exposure was calculated by weighted 
averaging over all orders not contaminated by telluric ab- 
sorption, and corrected for with the IRAfQ package RV . 
We found time and wavelength dependent variations in the 
blaze function between exposures, and cor rected for this us- 
ine: the method described in Winn et al. 1 (12004). We fitted 
the residual difference between an observed spectrum and 
an averaged template with an order 5 smoothed B-spline. In 
order to avoid removing the actual transit signal in the anal- 
ysis, we corrected the spectrum of order n with the averaged 
residual fits of orders n — 1 and n + 1. Finally, a template 
spectrum was created by a 2a clipped averaging of all blaze 
corrected object spectra. The Nal 'D' doublet region of the 
template spectrum is shown in Fig. [U 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation 



2.1.2 Subtraction of telluric features 

Telluric absorption lines originating from water and Na 
are dispersed unevenly within the stellar and interstellar 
Nal features (Fig[T|). Nightly variations in the strength of 
these telluric lines can overwhelm the planetary Nal absorp- 
tion signal, and must be removed in the spectral analysis 
process. Some transmission spectroscopy observations (e.g. 
Winn et al. 2004 : Narita et al. 2005; Redfield et al. 2008) re- 
moved their telluric contami nation by subtract i ng a smooth 
B sta r spectrum; others fe.g. lSnellen et al]|2008l : IWood et al.l 
l201lh performed the same subtraction using synthetic tel- 
luric spectra. We obtained a telluric template spectrum by 
averaging two observations of the B star HD 129116. For 
each WASP- 17 observation, the telluric spectrum, scaled 
and shifted by fitting for the telluric water dominated re- 
gions around the Nal 'D' doublet, was divided from the 
object spectrum. We note that the telluric scaling deviates 
from the airmass function at the start of the night, indicative 
of variations in the atmospheric water content. The telluric 
scaling correlates well with airmass for the second half of the 
night. We see no significant telluric Na lines in the telluric 
spectrum, and assum e that any telluric N a variations scale 
with water. Similarly, I Snellen et al.l (|2008l l pointed out that 
although telluric Na varies seasonally, it correlates well with 
water and oxygen over a single night. 



2.1.3 Flux extraction 

Following Equation 1 of I Snellen et al.l (|2008l l, we integrated 
the flux within a half-bandwidth of the line centroid, cor- 
rected against two adjacent reference continuum regions. 
We found that the blaze and reference continuum region 
corrections could not remove all systematic effects, but vari- 
ations in the interstellar Nal lines correlate strongly with 
that of the stellar Nal lines. The strong interstellar Nal ab- 
sorption doublet can be found 1 A (^ 50kms~^) red- ward 
of their stellar components (Fig.[T]). We can assume that the 
strength of the interstellar absorption features are invariant 
with time, as we do not expect the interstellar medium to 
be changing on the time scale of hours. The proximity of 
the interstellar to the stellar features makes them excellent 
references to correct for any additional, higher order system- 
atic variations. We hence extract the flux of both the stellar 
and interstellar absorption lines in our analysis. 

The Nal flux was calculated as the ratio between the 
sum of the contributing fluxes from each component of the 
doublet and that of the same regions in the template spec- 
trum. The line centroids of the stellar Nal 'D' doublet are 
located at 5888.89 (Di) and 5894.85 A (D2). The adjacent 
interstellar components are located at 5889.82 and 5895.80 

A. 

We find that the interstellar Di line is severely contam- 
inated by telluric water lines (Fig[T]). The line shows large 
time dependent variations, most likely due to the incom- 
plete subtraction of telluric lines blended within the stellar 
absorption feature. The relative level of telluric contamina- 
tion can be quantified by taking the sum equivalent width of 
the telluric lines within the stellar and interstellar regions; 
for the stellar Di and D2 lines the sum telluric equivalent 
widths are 0.0046 A and 0.0043 A, for the interstellar Di 
and D2 lines the sum telluric equivalent widths are 0.0092 A 
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Wavelength (A) 

Figure 1. Plotted are the Nal 'D' doublet regions of the template spectrum, with Di on the left panel, D2 on the right panel, 
as observed with Magellan. The broad stellar absorption features are resolved from the narrowe r interstellar components. Each 
line is also blended with additional interstellar and telluric Nal components (jWood et al.|[201ll V The minimum and maximum 
half-bandwidths of 0.60 and 0.90 A are filled in grey. The telluric absorption spectrum, as observed with a rapidly rotating B star, 
is over-plotted in blue. 



and 0.0008 A. The interstellar Di line contains more than 
twice the amount of telluric contamination as the stellar Di, 
D2 and the interstellar D2 lines. We therefore do not use the 
interstellar Di line in our subsequent analyses. 

For each exposure, we performed a Gaussian centroid 
fit to each line to remove any uncorrected wavelength shifts, 
then summed the flux around the fltted centroid position. 
The reference continuum regions to the left and right of 
the line centroid, for the stellar Nal Di line, are centred at 
5887 A and 5893 A respectively, and for the stellar and 
interstellar Nal D2 lines at 5893 and 5897 A respectively. 
These reference regions are not located directly adjacent of 
the stellar lines to avoid overlapping with the interstellar 
components. 

The use of the nearby interstellar Nal absorption lines 
as reference is facilitated by the large velocity difference 
between WASP- 17 and the sight line interstellar medium. 
The heliocentric mean velocity of WASP-17 is — 49.5kms~^ 
(fAnd erson et al. 2010), whilst we measure the sight line 
interstellar medium velocity to be — 0.5kms~^, consis- 
tent with the velocity of the local interstellar cloud (e.g. 
iRedfield k Linskvll2008l ). The separation between the stel- 
lar and interstellar Nal absorption features is ^ 1 A in wave- 
length space. To investigate the minimum stellar-interstellar 
feature separation required for this technique, we simulate 
a stellar, interstellar absorption pair with two equal height 
Gaussians of FWHM 0.4 and 0.2 A. Various separations be- 
tween the two Gaussians are tested, and the flux is extracted 
about the broader stellar Gaussian at 0.6 A bandwidth. We 
flnd that at separations less than 0.54 A, the influence of 
the interstellar line on the extracted stellar line flux becomes 
greater than our measurement error of 0.15%, corresponding 
to a minimum velocity separation requirement of 27.5 km s~^ 
between the target and the sight line interstellar medium. 

For a single exposure, the Nal flux, Frei, relative to the 
ratio between the stellar fluxes Fdi, stellar, i^D2, stellar and the 
interstellar flux Fd2, interstellar, is: 

Wl • FdI, stellar + W2 ' i^D2, stellar x 

i^rel ^ . (1) 

-t^D2, interstellar 



The weights wi and W2 are inversely proportional to the 
scatter of the lightcurves calculated for each stellar line in- 
dividually, according to 

i^Di ,stellar /^\ 

^Dl,rel — [^) 

-f^D2, interstellar 

and 

J-, -^D2, stellar /q\ 

^D2,rel — — • [o) 

-t^D 2, interstellar 

The flux of each individual component is calculated by 
their relative counts to the same regions in the template 
spectrum, corrected by the counts in the reference contin- 
uum regions. For example, the flux in the stellar Di line is 
derived from the counts in the object spectrum about the 
absorption line and reference region, Odi and Oref, and the 
counts in the same regions in the template spectrum, Tdi 
and Tref, by, 

17 _ ^Dl/^ref 

i'Dl, stellar — — T= • 

iDl/iref 

The error in the relative flux is derived from the average 
S/N within the bandwidth, divided by the square root of 
the number of spectral points included in the bandwidth. 



2.1.4 Modelling of in-transit line profile variations 

The strength of the absorption signal is strongly dependent 
on the bandwidth within which the flux is summed; the rela- 
tionship between the bandwidth and transit depth provides 
information on the properties of the planetary atmosphere. 
One limitation with our methodology of using the interstel- 
lar Nal lines as calibrators is that we must restrict the size 
of the bandwidth used. To avoid overlapping between the 
stellar and interstellar Nal components, the maximum half- 
bandwidth we employ is 0.90 A. The minimum bandwidth of 
0.60 A is determined by the line proflle variations in- transit 
that results from the distortions induced by the RM effect. 
The observed line proflle is a convolution between 



© 2012 RAS, MNRAS 000, [TH6] 



4 G. Zhou and D. D. R. Bayliss 



the intrinsic stellar absorption profile, a time varying ro- 
tational kernel that takes into account the occulted stel- 
lar disc, and the instrument response kernel. To model 
the distortions, the intrinsic Nal stellar absorption profile 
was approximated usi n^ a high-resolution solar spectrum 
(jDelbouille et al.lll990 V Solar rotational broadening is four 
times less than that of WASP- 17, whilst the instrumental 
broadening i n the solar template is 100 times less than our 
observations. iHirano et al.l (|2010l l provided a Gaussian ap- 
proximation to the analytical solution of the rotation ker- 
nel. Usin^ t heoretical V ban d limb darkening coefficients 
from lClaret'fc BloemenI (|201lh . with |Aridersone£aLl (|201lh 
stellar parameters for WASP- 17, we calculated a Gaussian 
rotational kernel with a standard deviation of 0.12 A. The 
kernel is modified at each time step to reflect the posi- 
tion of the planet over the stellar disc during transit. We 
modelled the MIKE instrument kernel with a Gaussian of 
0.052 A standard deviation, matching a spectral resolution 
of A/ A A = 48000 at 5890 A. The solar template is convolved 
at each time step, and the flux around the absorption line 
is summed and normalised against the out-of-transit flux. 

The upper panel of Fig. [2] shows the resulting transit 
lightcurve derived from the flux around Nal Di at various 
band widths. The distortions to the transit lightcurve, cal- 
culated as the absolute difference to the 2.0 A bandwidth 
lightcurve, become larger than our measurement error of 
0.15 per cent at bandwidths narrower than 0.60 A (Fig. [1 
lower panel). 

We therefore extracted the flux around the Nal 'D' dou- 
blet at seven half-bandwidths, from 0.60 to 0.90 A, evenly 
spaced at 0.05 A intervals. As an example, the 0.75 A half- 
bandwidth lightcurve is plotted in FigO and the fluxes are 
set out in Table [1] 



2.2 Lightcurve analysis 

The transit depth, D, is the ratio of the average out-of- 
transit flux to the average in-transit flux. Foot and Fint 
respectively ( Brown 2001 ), 



Foot 



(5) 



Exposures with mid-times set during ingress or egress are 
excluded. Overall, there were 20 out-of-transit points and 
11 in-transit points included in the calculations. The error 
in the transit depth, AD, is estimated from the standard 
deviation of the points in and out-of-transit, (Jint and aooT 
respectively. 



^ _ ( C"INt/ V A^INT 

y FiNT 



+ 



f aooT/VN( 
V Foot. 



OOT 



■ (6) 



where A^int and A'oot are the number of points averaged 
in and out-of-transit. The same analysis was performed 
at every bandwidth, the results of which are presented in 
Table O At 1.5 A bandwidth, we detect a transit depth 
of 0.58 di 0.13 per cent with 4.5cr confldence. We also per- 
formed the same lightcurve analysis on the fluxes extracted 
around the stellar Nal Di and D2 lines individually, ac- 
cording to Eq[2]and[3l The transit depths at 1.5 A band- 
width for the Di and D2 stellar lines are 0.53 =b 0.16 and 
0.66 ±0.10 per cent, respectively. The transit depths derived 
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Figure 2. The top panel plots the model transit lightcurves from 
fluxes extracted around the Nal Di line at bandwidths of 2.0 A 
(black solid), 0.7 A (green dashed), 0.5 A (red dot-dashed), and 
0.3 A (blue dotted). The bottom panel shows the absolute flux 
difference for the narrow bandwidths lightcurves, relative to a 2.0 
A bandwidth lightcurve. 
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Figure 3. An example lightcurve, extracted around the stellar 
Nal Di and D2 lines at 0.75 A half-bandwidth, is plotted. The 
red line plots the 0.58±0.13per cent transit depth fit. The darker 
shade marks the in-transit proportion of the lightcurve, whilst the 
lighter shade marks ingress and egress. 
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Table 1. Fluxes for the 0.75 A half-bandwidth lightcurve 





-f^rel 




CyomrQent 


0.544 


0.9918 


0.0096 


Out-of-transit 


0.554 


1.0073 


0.0083 


Out-of-transit 


0.564 


0.9998 


0.0078 


Out-of-transit 


0.573 


0.9987 


0.0081 


Out-of-transit 


0.583 


0.9992 


0.0077 


Out-of-transit 


0.592 


1.0021 


0.0071 


Out-of-transit 


0.602 


1.0021 


0.0068 


Ingress 


0.611 


0.9962 


0.0066 


Ingress 


0.630 


0.9941 


0.0066 


In-transit 


0.639 


0.9920 


0.0066 


In-transit 


0.649 


0.9935 


0.0073 


In-transit 


0.658 


0.9896 


0.0076 


In-transit 


0.668 


0.9871 


0.0069 


In-transit 


0.677 


0.9906 


0.0076 


In-transit 


0.686 


1.0011 


0.0078 


In-transit 


0.696 


0.9954 


0.0082 


In-transit 


0.706 


1.0026 


0.0079 


In-transit 


0.715 


0.9974 


0.0071 


In-transit 


0.741 


0.9929 


0.0084 


In-transit 


0.769 


0.9963 


0.0078 


Egress 


0.778 


0.9955 


0.0074 


Egress 


0.787 


1.0011 


0.0076 


Out-of-transit 


0.796 


1.0077 


0.0079 


Out-of-transit 


0.806 


0.9988 


0.0096 


Out-of-transit 


0.815 


0.9925 


0.0093 


Out-of-transit 


0.825 


1.0050 


0.0075 


Out-of-transit 


0.835 


0.9940 


0.0074 


Out-of-transit 


0.844 


1.0005 


0.0075 


Out-of-transit 


0.854 


0.9872 


0.0078 


Out-of-transit 


0.863 


0.9959 


0.0084 


Out-of-transit 


0.872 


1.0112 


0.0081 


Out-of-transit 


0.881 


1.0036 


0.0078 


Out-of-transit 


0.891 


1.0046 


0.0107 


Out-of-transit 


0.900 


1.0017 


0.0103 


Out-of-transit 


0.909 


0.9974 


0.0107 


Out-of-transit 



Table 2. The transit depths at various bandwidths. 



Bandwidth (A) 


Transit depth (per cent) 


1.2 


0.60 ±0.14 


1.3 


0.60 ±0.13 


1.4 


0.59 ±0.13 


1.5 


0.58 ±0.13 


1.6 


0.56 ±0.13 


1.7 


0.50 ±0.12 


1.8 


0.42 ±0.11 



from each individual NaT line is consistent with each other 
and with that derived from the weighted average lightcurve. 



3 DISCUSSION 

Our Magellan transmission spectroscopic measurements de- 
tected a 0.58 per cent transit, at 4.5cr confidence, in the NaT 
'D' doublet r efflon of WASP- 17 , which is consistent with 
the results of I Wood et alj (|201lh . This study confirms that 
WASP- 17b shows the strongest NaT absorption signal mea- 
sured to date. Given the low surface gravity and high equilib- 
rium temperature of the planet, WASP- 17b has the largest 




Bandwidth (A) 

Figure 4. Plotted are the various observed and modelled Nal 
transit depths. The Magellan detection presented in this paper 
is shown by red sha ded region. Blue diamonds plot the results 
of lWood et aP (l201lh. Gr een squares plot HD 209458b measure- 
ments from ISnellen etHI ( 2008) after scaling by t he WASP- 17b 
scale height. Yellow crosses plot the origina l ^Snellen e t al.l (|2008l ) 
measurements of HD 209458b. ThelBrownl ( 2001) model predic- 
tion, as calculated and scaled bv lWood et aimOll V is plotted as 
the dashed line. 

atmospheric scale height of all known transiting planets, fa- 
cilitating the large transmission spectroscopic signal. 

Interestingly, an extension of our analysis to narrower 
bandwidths could not replic ate the steep increase in absorp- 
tion strength as observed bv lWood et al.l(|201lh . It is unclear 
how the RM effect distortions to the absorption line profile 
were dealt with in this previous study, given that fluxes ex- 
tracted at 0.375 A half-bandwidth would have been severely 
affected. In addition, our results hint at a steep decrease in 
absorption strength for our broader bandwidths. This fea- 
tu re is potentially sim ilar to the steep decrease observed 
by IWood et alJ toU ) at 3.0 A bandwidth, which was at- 
tributed primarily to the presence of clouds in the upper 
planetary atmosphere. 

We can make a comparison between observations of HD 
209458b and that of WASP- 17b by scaling the former ac- 
cording to the differences in temperature and gravity be- 
tween the two planets. The transit depth is governed by the 
atmospheric scale height of the planet and the ra dius ra- 
tio between the planet {Rp) and the host star (R^) (|Brownl 
l200l|): 

_ 2ixRp{kT/^g) 
TvRi ' 

where T is the temperature, /a is the mean molecular weight 
of the atmosphere, and g is the surface gravity of the planet. 
Such scaling assumes that the atmospheric structure is con- 
sistent between the planets, and the change in contribution 
from secondary effects are negligible, accounting only for the 
difference in temperature and bulk density between the two 
planets. Taking a scale heig ht for HD 209458b of 500 km 
and /i of H2, we multiply the lSnellen et all (|2008l ) measure- 
ments by a factor of 2.4, the resulting estimates are plotted 
in Fig. m Our WASP- 17b measurements match the scaled 
HD 209458b atmosphere remarkably well. 

Similarly, the transit depths derived and scaled from 
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models of 'Brown (200l") bv lWood et aD ('2011') are also plot- 
ted. The model prediction lies well above our measured sig- 
nal. All transmission spectroscopy of transiting planets to 
date have consistently shown that the Nal signal is weaker 
than that of a cloudless model in local thermodynamic equi- 
librium. The close-in orbits of these 'hot-Jupiters' mean that 
their atmospheres are highly irradiated, causing depletion of 
Na to form Na^ . Photoionisation occurs in the upper levels 
of the atmosphere, whilst the line wings are formed in the 
denser lower atmosphere. The result of neutral Na deple- 
tion is ther efore a reduction in the depth o f the line cores 
fe.g. Barman et al.ll200i iFortnev et al.ll2003i ). The presence 
of silicate and iron condensates in the upper atmosphere 
(jSudarskv et al.ll2QQ0l . l2003l : iFortnev et al.|[2003l ). combined 
with the slanted geom etry inherent to transmission spec- 
troscopy (|Fortnevli2005l ). also act to decrease the absorption 
signal. 

We note that this study presents a ground based spec- 
troscopic Nal detection with the smallest telescope to date, 
in a S/N regime significantly lower than previous observa- 
tions. This shows that the use of interstellar absorption fea- 
tures as reference is a robust technique that removes the 
majority of systematic effects. iLangland-Shula et al.l (|2009l l 
noted that better telluric subtractions could be achieved by 
taking interleaved smooth spectrum star observations, re- 
ducing the reliance on scaling a single telluric template to 
fit the observed spectrum. WASP- 17 is by far the faintest 
target pursued by the transmission spectroscopy technique, 
the positive detection opens up the possibility to explore 
other similarly inflated targets. 
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